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ABSTRACT 

The present acceleration of the Universe leads to the formation of a cosmological 
future event horizon. We explore the effects of the event horizon on cosmological 
backreaction due to inhomogeneities in the universe. Beginning from the onset of the 
present accelerating era, we show that backreaction in presence of the event horizon 
causes acceleration to slow down in the subsequent evolution. Transition to another 
decelerating era could ensue eventually at a future epoch, ensuring avoidance of a big 
rip. 



1 INTRODUCTION 



There exists overwhelming observationa l evidence for the 
present accele r ation of the Universe (IPerlmutter 1 19981 ; 
iKowalskil [20081 ; iHickenl 120091 ; ISeikel fc Schwarzll2009h . The 
accelerating universe leads to a future event horizon from 
beyond which it is not possible for any signal to reach 
us. On the other hand, observations also tell us that 
our Universe is inhomogeneous up to the scales of su- 
per clusters of galaxies. The idea that backreaction orig- 
inating from density inhomogeneities could lead to modi- 
fications in evolution of the universe as described by the 
background Friedmann- Robertson- Walker (FRW) metric at 
large scales has gained popular i ty in recent years (Buchert 
200d; IBuchert fc Carfora~ll200Sl; Iwiegand fc Buchert l |20ld; 



Zalale tdinovlll992l ; Buchert fc Carfora ll2003l;lRasaner]20o4 
| 200Sl.|20ld;)Kolb et al.ll2005l.l200Sj ; |Paraniape fc Singhlboogt 
Singbl l201l(; iMattson fc MattsonI [20ld : IWiltshirej l2007al lbl; 
IGasperini et al.l 120090 . Here we show that backreaction in 
the presence of the cosmological event horizon could have a 
remarkable consequence of ushering in another decelerated 
era beyond the present accelerating epoch. 

In spite of nume rous creative ideas proposed for th e 
present acceleration |Sahnil |2004| ; ICopeland etafl l200rj ). 
there is still a lack of convincing explanation of this phe- 
nomenon. The simplest possible explanation provided by a 
cosmological con stant is endowed with conceptual problems 
l|Weinberd[l989l l. Alternative mechanisms based on either 
modifications of the gravitational theory, or invoking extra 
fields suffer from the coincidence problem, as to why the 
era of acceleration begins around the same era when the 
Universe becomes structured. The ultimate fate of our Uni- 
verse remains clouded in considerable mystery. Backreaction 
from inhomogeneities provides an interesting platform for 
investigating this issue without invoking additional physics, 
since the effects of backreaction gain stength as the inhomo- 
geneities develop into structures around the present era. 

Approaches have been developed to calculate the ef- 
fect of inhomegeneous matter distibution on the evolu- 



tion of the Universe ijBuchertJ |200(]I : IZalaletdinovl 1 19921 : 
IKolb et all 120051 ). Arguments in favou r of the viabilit y 
of backreaction seem rather compelling |Kolb et al.l 12008). 
though there exists debate on the impact of inhomogeneities 
on observables of an ov e rall homogeneous FRW model 
JParaniape fc Singhl 120081 ; llshibashi fc Waldl l200rj ; ISinghl 
[20110 , and on the magnitude of backreaction modulated by 
the effect of shear between overdense and underdense regions 
|Mattson fc Mattson l201(ih. Using the framework f omulated 
by Bucher t jBuchertf EoOO'l it has been shown ()Rasanenl 
|2004 |200S|, l20irj) that backreaction could lead to an accel- 
erated expansion during the present epoch. A notable appli- 
cation of the formalism has been developed by Wiltshire 
showing an apparent volume acceleration of the universe 
based on the different la pse of time betw een the underdense 
and overdense regions ((Wiltshire! |2007al fb|) . Further, gauge 
invariant averages in the Buchert framewor k have also been 
constructed recently |Gasperini et al.ll2009T ). 

While upcoming observations may ultimately de- 
cide whether backreaction from density inhomo geneities 
drives the present acceleration, the above studies (IBuchertj 
200d; IBuchert fc Carfora 1 1 20081: IWiegand fc Buchert | | 201 
Buchert fc Carfora I 12003; iRasaneni 12004 |200S|. I20T 



Kolb et all 120051 . 120081 ; IWiltshirej l2007al lbl: IGasperini et all 
2009) have highlighted that backreaction could be an im- 



portant ingredient of the evolution of our Universe. Here we 
explore this issue with a fresh perspective, viz., the impact 
of the event horizon on cosmological backreaction. The cur- 
rently accelerating epoch dictates the existence of an event 
horizon since the transition from the previously matter dom- 
inated decelerating expansion. Since backreaction is evalu- 
ated from the global distribution of matter inhomogeneities, 
the event horizon demarcates the spatial regions which are 
causally connected to us and hence impact the evolution 
of our part of the Universe. In the present work we inves- 
tigate the consequences of backreaction in presence of the 
horizon. Such an approach has remained unexplored in pre- 
vious studies of backreaction. It may be noted that the for- 
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malism of back reaction llBuchertll2000l ; [Buchert fe Carforal Q-p is given by 



l2008l ; IWiegand fc Buchert Il2010l ) has been criticized on the 
grounds that the average is taken on a space like hypersur- 
face, whil e observations are made along and inside the past 
light cone llshibashi fc Waldl ((2006). Our present analysis, by 
considering an effect due to the event horizon, introduces an 
element of light cone physics from a somewhat different per- 
spective. We show that backreaction with the event horizon 
could lead to the possibility of transition to a decelerated 
future era. 



2 THE BACKREACTION FRAMEWORK 

In the framework develo p ed by Buchert (Buchert |2000| ; 



iBuchert fc Carfora I 120081 ; IWiegand fc Buchert I I20T0I ) for 



the Universe filled with an irrotational fluid of dust the 
spacetime is foliated into flow-orthogonal hypersurfaces fea- 
turing the line-element ds 2 = — dt 2 + gijdX x dX\ where 
the proper time t labels the hypersurfaces and X r are 
Gaussian normal coordinates (locating free-falling fluid el- 
ements or generalized fundamental observers) in the hyper- 
surfaces, and g lJ is the full inhomogeneous three metric 
of the hypersurfaces of constant proper time. For a com- 
pact spatial domain T> whose volume is given by \T>\ g — 
f^d^g where dp, g = y/^gjt, A 1 , X 2 , X 3 )dX 1 dX 2 dX 3 , 

/ \T>\ \ 1/3 

the scale factor ap(t) — ( p | j encodes the average 
stretch of all direc tions of the domain. The Einstein equa- 
tions then lead to |Bucherdl2000l : IBuchert fc Carfora llioosl : 
IWiegand fc Buchert IteoTcf ) 



3^ 

ml 



-4%G{p) TI + Qt> + A 
&xG(p) v -±(1l) v -±Qv+A 



(1) 



= d t (p) v + 3Hp (p) 7 



where the average of the scalar quantities on the domain T> is 

defined as (/) v (t) = ^ ^Xf^ = fv f d *>, 

and where p, 1Z and Hp denote the local matter density, the 
Ricci-scalar of the three-metric gij, and the domain depen- 
dent Hubble rate Hp = dp jap respectively. The kinemati- 
cal backreaction Qp is defined as Qp = | ((^ 2 ) IJ — (6)p) ~~ 
2cr|>, where 6 is the local expansion rate and a 2 = XjlaijO 11 
is the squared rate of shear. Qp encodes the departure from 
homogeneity. 

The "global" domain T> is assumed to be separated into 
subregions Tt which themselves consist of elementary space 
entities T^ that may be associated with some averag- 
ing length scale, i.e., T> — UeTe, where Te = U a Tf and 



T^ n t; 



for all q 7^ j8 and I 7^ m. Analogous to the 
scale factor for the global domain, a scale factor ae for each of 
the subregions Te can be defined such that \D\ g = Y^e \J~e\ g , 
and hence cip = ^2 e ^ e i a £> where \e i — l-T^ j s /|2?ij 9 is the 
initial volume fraction of the subregion Te- The average of 
the scalar valued function / on the domain T>, may then 
be split into the averages of / on the subregions Te in the 
form, (f) v = El^ls X E /tW fdUg = £-M/)jv where 

l a £ I 

A« = \J r e\g/\1- > \g, is the volume fraction of the subregion Te- 
Due to the (6) v term, the expression for the backreaction 



Qp = ^^Af + AfA m (He — H„ 



(2) 



where, Qe and He are defined in Te in the same way as 
Qp and H-p are defined in D. The shear part (cr 2 ) T is 
completely absorbed in Qe whereas the variance of the lo- 
cal expansion rates (O 2 )^ — {6)%> is partly contained in Qe 
but also generates the extra term 3 Ef^m Af A m (He — H m ) 2 . 
This is because the part of the variance that is present in 
Qe, namely (0 2 ) T — (8)% t only takes into account points in- 
side Te- To restore the variance that comes from combining 
points of Te with others in T m , the extra term containing 
the averaged Hubble rate emerges. Note here that the above 
formulation of the backreaction holds in the case when there 
is no interaction between the overdense and the underdense 
subregions. 

Now from Eq.([TJ one gets 



dp _ aejt) 
a-D ae(t) 



Following 



the 



+ ^XtXm (Hi - H m f (3) 

of 



sim plifyin g assumption 
Ref. (|Wiegand fc Buchertl |201Ch . (which captures the 
essential physics) we work with only two subregions. Club- 
bing those parts of T> which consist of initial overdensity 
as M (called "wall"), and those with initial underdensity 
as £ (called "void"), such that T> = M U £, one obtains 
Ht> = \mHm + XeHs, with similar expressions for (p)^, 
and {TZ)^, and 



dp 



\m—— + A^— + 2\m^e(Hm 



ajvi 



(I E 



He) 



(4) 



Here \ M + A £ = 1, with X M = \M\/\T>\ and A £ = |£|/|2?|. 
Since the global domain T> is large enough for a scale 
of homogeneity to be associated with it, one can write 
\V\ g — J v \f—g d [i X w f(r)a%(t), where f(r) is a function 
of the FRW comoving radial coordinate r. It then follows 

that ax> ~ (^>p) aF ' an< i hence, the volume average 
scale factor ax> and the FRW scale factor of are related by 
a-D ~ cfclf, where cf is constant in time. Thus, Hf ~ H-p, 
where Hp is the FRW Hubble parameter associated with 
T>. Th ough in general Hp and Hf could differ on even large 
scales (jWiegand fc Buchert l201oh . the above approximation 
is valid for small metric perturbations. 



3 EFFECT OF EVENT HORIZON 

We now come to the central issue of the paper, as to what 
happens to the evolution of the universe once the present 
stage of acceleration sets in. Note henceforth, we do not need 
to nece ssarily assume that the acce l eration is due t o backre- 
action l|Wiegand fc Buchert |[201ol : |R asane nl 120041 ). For the 
purpose of our present analysis, it suffices t o consider the ob- 
serve d accelerated phase of the universe |Seikel fc Schwarzl 
l2009h that could occur due to any o f a variety of mechanisms 
(|Sahnill2004l ; ICopeland et~"ai1l2006h . Given that we are un- 
dergoing a stage of acceleration since transition from an era 
of structure formation, our aim here is to explore the sub- 
sequent evolution of the Universe due to the effects of back- 
reaction in presence of the cosmic event horizon. Though 
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Figure 1. The dimensionless global acceleration parameter 
is plotted versus time (s) assuming a constant horizon. The 

values for the various parameters used are (i) a = 0.995, = 0.5, 
(ii) a = 0.999, = 0.6, (iii) a = 1.0, /3 = 0.5, and (iv) 
a = 1.02, = 0.66. 



in general, spatial and light cone distances and correspond- 
ing accelerations could be diff erent, as shown explicitly i n 
the framework of LTB models l|Boleiko fc An dcrsson 2008), 
an approximation for the event horizon which forms at the 
onset of acceleration could be defined by 



dt' 

Ou(t') 



(5) 



in the same spirit as ap = Ci?ap. 

Following the Buche rt framework |Buchertl |2000| : 
IWieeand fc Buchert I I20T0I ) as discussed above, the global 
domain T> is divided into a collection of overdense regions 
M — VJjM\ with total volume \M\ g = Ylj l-M 3 |gj an d 
underdense regions £ = Uj£ 3 with corresponding volume 
\£\ g = J^. \£ 3 \ g - Assuming that the scale factors of the re- 
gions £ J and M 3 are respectively given by aej = eg, t a and 
o-Mj = CMjt 13 , where a, j3, ce, and CMj are constants, one 
has 



3 3 ,3a 

a-e = c £ t ■ 



3 _ 3 ,3/3 
a M — c M l 



(6) 



where cf 



\£i\ 



is a constant, and similarly for 



cm. The volume fraction of the subdomain M is given by 



\T>\ 9 



which can be rewritten in terms of the corre- 



sponding scale factors as \m = 



Since an event 



horizon forms, only those regions of T> that are within the 
event horion are accessible to us. Hence, in this case an ap- 
parent volume fraction \m h given by \m h = ^gir-r^ is 
introduced. From eq.© it follows that 



(7) 



where c Mh = 3c^|A1i|g/47r is a constant. Normalizing the 
total accessible volume in the presence of the event horizon, 
we can write 

\s h = 1 - \M h (8) 
where \e h is the apparent volume fraction for the subdomain 



£ . It hence follows that the global acceleration equation ([4]) 
is now given by 



ap 



1) 



t 2 



„3 ,3/3 
C M h t 



a(ct — 1) 



+2 



r 3 i 3/3 



C M h t 



3/3 



(9) 



In order to obtain the future evolution of the universe 
with backreaction in presence of the event horizon, one has 
to solve the above equation for the scale factor with the 
event horizon given by Eq.Q. In what follows we will 
eventually obtain numerical solutions of the above integro- 
differential equations. However, it is first instructive to ob- 
tain some physical insight of the evolution by taking recourse 
to a simple approximation. 

To this end, let us for the moment model the onset of 
the present acceleration of the Universe by an exponential 
expansion, keeping our analysis close to observations. Specif- 
ically, we set ap oc e Hl3t in Eq.Q only. (We will see later 
that this rather crude approximation does indeed give rise 
to results that are qualitatively similar to the ones obtained 
through numerical analysis). Using Hf = Hp, where Hf 
is the FRW Hubble parameter associated with T>, it follows 
that rh = Hp 1 , a constant which we substitute in Eq.©. 
With this substitution, the global acceleration 'dp vanishes 
at times given by 



3/3 



4(/3- 



(3/3- 



1) ± d (3/3 - a - I) 2 +8a(a- 1) 



factor of the "wall" grows as v , where 



(10) 
1/2 < 



Eq. (|10p corresponds to real time solutions for 



The scale 
P < 2/3. 

a >§[(/} + 1) + 2^20(1-/3)] . 

Now, let us consider the following two cases separately: 
Case I: a < 1 and /3 ^ 2/3. There exist two real solutions 
(| 10f) corresponding to two values of time when the global ac- 
celeration vanishes. In Fig.l we plot a dimensionless global 
acceleration parameter °° » with time using eqiH). The 

curves (i) and (ii) correspond to this case showing that the 
Universe first enters the epoch of acceleration due to backre- 
action, which subsequently slows down and finally vanishes 
at the onset of another decelerating era in the future. Case 
II: ot^tl and /3 < 2/3. From {TOjl it follows that there is only 
one real solution (minus sign for the square root). This case 
models the Universe which accelerates due to some other 
mechanism (not backreaction), but subsequently enters an 
epoch of deceleration due to backreaction of inhomogeneities 
in the presence of the event horizon (see curves (iii) and (iv) 
of Fig.l). 

The plots in Fig.l have been done taking the standard 
values of the parameters rh = H^, 1 — 4.36 x 10 17 s while 
choosing the appropriate range for the parameters a and /3, 
as given in the figu re caption. Based on the N- body simula- 
tion values used in l|Wiegand fc Buchert |[201ol ) we also take 
Aa-i^o = 0.09. Using the relation zt = exp [Hp (to — tr)] — 
1, where tr corresponds to the transition time in the past, 
the red-shifts for the transition could be estimated. For ex- 
ample, for the data used in curve (i), the transition from 
deceleration to acceleration occurs at zt — 0.844, and for 
curve (ii) we have zt — 0.914 (which are close to the ACDM 
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le+017 2c+017 3e+017 4e+017 5e+017 6e+017 
t 

Figure 2. The dimensionless global acceleration parameter 
- ai ^'2 is plotted versus time (s) as obtained through numerical 
integration, with the 'initial condition of go = —0.7. The values 
for the various parameters used are chosen to be the same as in 
the corresponding plots of Fig.l. 



value for the standard transition redshift (Mel chiorri et all 
[2003)). 

We now study the acceleration equation ^ numerically 
without assuming a priori any behaviour for the horizon. 
Keeping with the spirit of our analysis, we assume that the 
Universe has entered the accelerated stage and thus a cosmic 
event horizon has formed. This ensures that ru defined by 
((5| will be finite valued, enabling us to replace the integral 
equation ((SJ) by 



? 

an 



(11) 



Thus, the evolution of the scale factor is now governed by 
the set of coupled differential equations @ and We 
numerically integrate these equations by using as an "initial 
condition" the observational constraint qo w —0.7, where go 
is the current value of the deceleration parameter, and using 
the solution for the scale factor plot the global acceleration 
versus time in Fig. 2 (thus all the curves in Fig. 2 are set to 
intersect at the point (to, go))- The values of the other pa- 
rameters including a and /3 are chosen to be the same as 
in the corresponding curves of the exponential case. As can 
be seen from Fig. 2, the nature of the plots is quite similar 
to the ones that were obtained in the case assuming a con- 
stant event horizon, with the a > 1 curves signifying only 
one transition between acceleration and deceleration in the 
future. The differences in the various slopes and also in the 
scale for the dimensionless global acceleration parameter in 
the two cases arise as a result of the approximation of con- 
stant horizon used in the former, as well as due to the choice 
of the condition go ~ —0.7 used in the latter. 



4 CONCLUSIONS 

To summarize, in this work we have explored the effect 
of backreaction due to inhomogeneities on the evolution 
of the Universe undergoing present acceleration. We have 
shown that the presence of the cosmic event horizon causes 
the acceleration to slow down significantly with time. Our 
results indicate the fascinating possibility of backreaction 



being responsible not only for the present acceleration as 
show n in earlier works ( Wie gand fc Buchert Il20l"ol ; iRasanen] 
2010), but also leading to a transition to another deceler- 
ated era in the future. Another possibility following from 
our analysis is of the U niverse currently accelerating due t o 



but with backreaction dBuchert 


2000; Buchert & Carfora 


I200SL IWiegand & Buchert 201fJ 


) later 


causing acceleration 



to slow down. Our prediction of the future slowing down 
of acceleration seems to fit smoothly with the earlier era of 
structure formation and the transition to acceleration in the 
standard ACDM model, as shown here (transition red-shift 
z T w 0.8). 

Before concluding, in context of the formalism used in 
the present work it may be worthwhile to recapitulate some 
of the present debate in t he literature regarding averaging on 
a space like hypersurface dBuc hert 2000; [Buchert fc Carforal 
2008; Wiegand fc Buchertll2010r ) as compared to taking the 
average on the past light cone. The usefulness of the ex- 
pansion rate averaged on any hypersurface is determined by 
relating it to observed quantities. It has been observed that 
the redshift and distance can be expressed in terms of the 
average geometry alone, p rovided that th e contribution of 
the null shear is neglig ible (|Rasanenll2009h. Observationally, 



the shear is known to be indeed small ( Munshi et al.|[20Q. 



Nonetheless, it has been claimed that neither averaging on a 
constant time hypersurface nor light cone averaging is easy 
to connect with the obser vations correspondin g to parame- 
ters of the ACDM model (|Kolb fc Lambll2009h . The task of 
developing a procedure for light cone averaging is an ambi- 
tious program and till date there is no standard formalism 
to do so. In a recent paper three different types of ligh t 
cone averaging have been proposed (|Gasperini et al.ll201lf ). 
though much work remains to be done in order to apply 
their technique to the problem of cosmic acceleration. On 
the other hand, our present work introduces an element of 
light cone physics from another perspective by considering 
an effect due to the event horizon. 

Finally, it may be noted that though the event hori- 
zon is observer dependent, it follows from the symmetry of 
the equations Q and {5J| that our analysis leads to simi- 
lar conclusions for a 'Void" centric observer, as it does for a 
"wall" centric one. Note also that our analysis is valid while 
the event horizon exists. Hence, if the acceleration vanishes 
at some epoch in the future, one needs to consider backre- 
action without the event horizon beyond that epoch. The 
scales which crossed ouside the horizon earlier, will begin 
re-entering with backreaction from their associated inhomo- 
geneities impacting the evolution. Such a scenario is some- 
what reminiscent of the horizon crossing of modes during 
inflation in the early universe, and their subsequent reentry 
with rich cosmological consequences. In the present context, 
the impact on the global evolution of the reentering scales 
needs to be studied further. Moreover, it would be worth- 
while to investigate effects of the cosmic event horizon on 
other models of backreaction in order to make more generic 
predictions of observational interest. 
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